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Abstract 



We discuss the profile of neutrino masses and mixings in models with large 
extra dimensions when right handed neutrinos are present in the branes along 
with the usual standard model particles. In these models, string scale must 
be bigger than 10 8 GeV to have desired properties for the neutrinos at low 
energies. The lightest neutrino mass is zero and there is oscillations to sterile 
neutrinos that are different from other models with the bulk neutrino. 



I. INTRODUCTION 

Last year has seen an explosion of interest and activity in theories with large extra di- 
mensions [|1]-§J. It has been realized that extra dimensions almost as large as a millimeter 
could apparently be hidden from many extremely precise measurements that exist in particle 
physics. What makes such an idea exciting is the hope that the concept of hidden space 
dimensions can be probed by collider as well as other experiments in not too distant future. 
Furthermore, it brings the string scale closer to a TeV in some scenarios, making details of 
string physics within reach of experiments. On the theoretical side, recent developments in 
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strongly coupled string theories have given a certain amount of credibility to such specu- 
lations in that the size of extra dimensions can be proportional to the string coupling and 
therefore in the context of strongly coupled strings such large dimensions are quite plausible. 

This new class of theories have a sharp distinction from the conventional grand unified 
theories as well as old weakly coupled string models in that in the latter case, most scales 
other than the weak scale and the QCD scale were assumed to be in the range of 10 13 to 
10 16 GeV. That made it easier to understand observations such as small neutrino masses 
and a highly stable proton etc. Now that most scales are allowed to be small in the new 
models, the two particularly urgent questions that need to be answered are why is the proton 
stable and why are neutrinos so light. It has been speculated that proton stability may be 
understood by conjecturing the existence of U(l) symmetries that forbid the process. On 
the other hand no such simple argument for understanding small neutrino masses seems to 
exist. The familiar seesaw mechanism 1(J is not implementable in simple versions of these 



models where all scales are assumed to be low. So understanding lightness of neutrinos is 
therefore a challenge to such models. 

One approach to this problem discussed in recent literature ]7||| is to use neutrinos 
that live in the bulk (they are therefore necessarily singlet or sterile with respect to normal 
weak interactions) and to observe that their coupling to the known neutrinos in the brane 
is inversely proportional to the square-root of the bulk volume. The neutrino mass (which 
is now a Dirac mass) can be shown to be ~ h ^ , where M* is the string scale. If the 
string scale is in the TeV range, this leads to a neutrino mass of order 10~ 4 eV or so. This 
therefore explains why the neutrino masses are small. A key requirement for small neutrino 
masses is therefore that the string scale be in the few TeV rangeQ. Furthermore, these 
models also generally predict oscillations between the modes of the bulk neutrino and the 
known neutrinos. Thus if one wants an explanation of the solar neutrino deficit via neutrino 
oscillations, this implies that there must be at least one extra dimension with size in the 
micrometer range IHJ. The latter provides an interesting connection between neutrino 
physics, gravity experiments searching for deviations from Newton's law at sub-millimeter 
distances as well as possible collider search for TeV string excitations. It then follows that 
larger values of the string scale would jeopardize this simple explanation of the small neutrino 
mass. So the extent that larger values of string scale are also equally plausible as the TeV 
value, one might search for alternative ways to understand the small neutrino masses. 

It is the goal of this paper to outline such a scenario and study its consequences. The 
particular example we consider illustrates this scenario in models with large extra dimensions 
and generic brane-bulk picture for particles, where the string scale is necessarily bigger 
(> 10 8 GeV or so) and solar or atmospheric neutrino oscillations require at least one extra 
dimension be in the micrometer range. The new ingredient of the class of models we discuss 
here is that we include the right-handed neutrino in the brane and consider the gauge 
interactions to be described by a left-right symmetric model. In these models, the left- 
handed neutrinos are not allowed to form mass terms with the bulk neutrino due to extra 
gauge symmetries; instead it is only the right-handed neutrino which is allowed to form 
mass terms with the bulk neutrinos. This leads to a different profile for the neutrino masses 
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and mixings. In particular, we find that in this model, the left-handed neutrino is excatly 
massless whereas the bulk sterile neutrinos have masses related to the size of the extra 
dimensions and it will be of order ~ 10~ 3 eV if there is at least one large extra dimension 
with size in the micrometer range. A key distinguishing feature of this model from the 
existing ones is that the string scale is now necessarily much larger than a TeV. We also find 
that the pattern of the neutrino oscillations is different from the previous case. 

As mentioned, the minimal gauge model where our scheme is realized is the left-right 
symmetric models where the right handed symmetry is broken by the doublet Higgs bosons 
Xr(1,2, 1). The notation we follow is that the three numbers inside the parenthesis cor- 
respond to the quantum numbers under SU{2) L x SU{2) R x U(1)b-l- We do not need 
supersymmetry for our discussion and will therefore work within the context of nonsuper- 
symmetric left-right models. 

To set the stage for our discussion, let us start with a review of the neutrino mass 
mechanism in models with large extra dimensions discussed in Ref. J7J. The basic idea is 
to include the coupling of bulk neutrino ^(x^, (which is a standard model singlet) to 
the standard model lepton doublet y — 0). The Lagrangian that is responsible for the 

neutrino masses in this model is: 

S = k J d 4 xLHu B (x,y = 0) + J d 4 xdyi , B(x,y)T 5 d 5 i' B (x,y) + h.c. (1) 
Writing the four component spinor v B = ( . Pb 2 * j , we get for the neutrino mass matrix: 

<ww(7f ){z) p) 

This is a compact way of writing the KK excitations along the fifth dimension. Notice that 
the a/2 in the off-diagonal term appears to compensate the different normalization of the 
zero mode in the Fourier expansion of the bulk field in terms of smny/R and cos ny/R (R 
being the radius of the fifth dimension). Also, notice that only the last terms may couple 
to the pure four dimensional fields, and we always may choose these modes to have positive 
KK masses. In the above equation k embodies the features of the string scale and the radius 
of the extra dimension i.e. k ~ jr~- F° r KV ^ A*o where /io = -R^ 1 , we get the mixing of 
the v e with the bulk modes to be 



>eB 



dd — k 2 v 2 



(3) 



Substituting the eigenvalues of the operator <9 5 = n/i , we get for the n-th KK excitation a 
mixing tan6 n ~ ^ where £ w :y ^p. This expression is same as in [|7|||. 

Important point to note here is that since v e has a mass of kv, present neutrino mass 
limits lead to an upper limit on k and hence on the string scale. For instance, if we choose 
the present tritium decay bounds |TT| of m Ue < 2.5 eV, we get M* < 10 7 GeV. This bound 



gets considerably strengthened if we further require that the solar neutrino puzzle be solved 
via v e — v s oscillation (where we have called the typical excited models of the bulk neutrino 
as "sterile neutrinos"). The reason is that this requires Am 2 ~ 10~ 5 eV 2 and in the absence 
of any unnatural fine tuning, we will have to assume that that m Ue ~ m Us ~ 1CT 3 eV. This 
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implies that M* < 10 TeV and the bulk radius is given by R^ 1 ~ 1CT 3 eV implying R ~ 0.2 
mm. Another prediction of this model is that the mixing of the v e with the bulk neutrinos 
goes down like (1/n) where n denotes the level of Kaluza-Klein excitation. 

Let us now proceed to the new case we are considering. The part of the action relevant 
to our discussion is given by: 

S = J d A x[KLxL^B{.x J y = 0) + KRxR v B{x,y = 0) + hL<j)R} + J d^xdyVB^ d^vs + h.c. (4) 

where L T = (z/ e£ , e£) and R T = [y eR , e R ) and is the bidoublet Higgs field that breaks the 
elctroweak symmetry and gives mass to the charged fermions. We assume that the SU(2) R 
gauge group is broken by < Xr >= V R with < x°l > = 00 an d the SU(2)l x U(1)y is broken 
by Diag < <p >= (v,v'). The profile of the neutrino mixing matrix in this case is given by 

ipeL V*BL V' BL ) | KV R ] ( ^J R ) (5) 

where hv is the usual Dirac mass term present in the models with the seesaw mechanism 
and normally assumed to be of order of typical charged fermion masses (we will also make 
this plausible assumption that hv is of the order of a few MeV's for the first generation 
which will be the focus of this article). 

We can now proceed to find the eigenstates and neutrino mixings. First point to note is 
that this matrix being 2x3 has one zero eigenvalue corresponding to the state [|T2|] 

u = (cosC v eL - sinC vqbl) (6) 

where sin( = , hv . If we want the lightest eigenstate to be predominantly the 

y/{KV R ) 2 + (hv) 2 

electron neutrino so that observed universality in charged current weak interaction is main- 
tained, we must demand that kv r ^> hv. Since k ~ , this constraint will imply a 
constraint on the string scale M* . 

Let us see under what circumstances this condition is satisfied. Since the Dirac mass 
hv ~ few MeV's, we would like kv r ^> few MeV's. Let us assume that there is one extra 
dimension with large size (of order milli- meter and denoted by R\) and all other extra 
dimensions have very small sizes, assumed to be equal. Then the observed strength of 
gravitational interaction implies the relation 

M 2 pe = M* n+2 R n ~ 1 R 1 (7) 



2 In general, in the left-right model, there is a coupling of the form Mqxl4>XR which induces a vev 
for the x\ °f or der ' wn i cn f° r a choice of Mq = v w k gives a < \\ >~ 0.1 MeV. We 

assume Mq to be smaller so that this contribution to the neutrino masses is negligible. It could 
also be that the symmetries of the charged fermion sector are such that they either forbid such a 
coupling or give such a vev to (j> that this term does not effect the value of the potential at the 
minimum. We thank the referee for pointing this out. 
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where R\ is the largest dimension and the rest of the R's are small as required by neutrino 
physics. For simplicity let us identify the right handed symmetry breaking scale, the string 
scale and the inverse radii of the small dimensions R~ l . Then we have the approximate 
relation that 



1 



M* 



K 



^M* n R n - 1 R 1 M Pl 



(8) 



kv r > few MeV (say 100 MeV), implies that M* ~ v R « R- 1 > 10 8 5 GeV. In fact it is not 
hard to see that to satisfy the relation in Eq. (]?]), the radii of the "small" compact dimensions 
must be also of order M* _1 and that of the large dimension is of course in the sub-millimeter 
range. While this is a generic possibility, one can of course make many variations on this 
general theme. The results of this paper are not effected by such variations. Thus it appears 
that our scheme prefers a high string scale in contrast with the earlier proposal 0. 

We can now look at the rest of the neutrino spectrum arising from the KK excitations 
of the bulk mode as well as the right handed neutrino. They can be studied by looking at 
the "2 x 2" matrix after extracting the zero mode discussed above. Defining the orthogonal 
combination to z/ as z> l, we have: 



f ~ ./ J kv r 




(9) 



where z> l = cos( vobl + sm( v e L and we have used the approximation that kv r ^> hv. In 
terms of this short handed notation, the characteristic equation of the Dirac mass matrix 
above is easily computed to be 



m 



dl) 



m 



2 2 i 
KV R + 



2m 2 K 2 v R 



0. 



(10) 



This, in a manner similar that noted by Dienes et al. [|7| translates into the transcendental 
equation 



7tk v R R cot(7rm n R) 



(11) 



where m n is the mass of the n-th KK state. The tower of eigenstates, symbolically denoted 
by v n L, are exactly given by 



1 

Vn 



vol + 



with the normalization factor, r] n , given as 



V2m 2 n 
m 2 — d 2 



2miR 4 



BL 



n 2 = i + y 



(12) 



(13) 



As long as kv r R 3> 1, all the mass eigenvalues, m n upto m n ~ kv r satisfy cot{nm n R) ~ 
0. Therefore, the masses for all these states are ~ (2n — 1)/2R. Thus the effect of the mixing 
of the brane righthanded neutrinos is to shift the bulk neutrino levels below and near kv r 
downward. This is similar to what was noted in the case without the v R in the brane @ 
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for all states. On the other hand, using the same equation (Eq. (11)) it is easy to see that 
the states much heavier than kvr have masses ~ n/R, which also means that they basically 
decouple from the light sector. For the eigenstates in the middle, those with masses just 
beyond kvr, their mixing with the lightest states strongly suppress their contribution to 
z> L, by an amount ~ 1/kvrR, as it could be checked from Eq. ([13]), by summing over those 
elements for which k <C kvrR. 

There are, of course, similar results for the right handed states which involve the right 
handed bulk neutrinos. In fact, they are identified as the negative solutions to the charac- 
teristic equation, since the same mass matrix gives the masses of both left and right handed 
sectors and those are degenerate. 

For discussing neutrino oscillation, the left handed eigenstates are the only ones relevant. 
For this purpose let us write down the weak eigenstate v e in terms of the mass eigenstates 
through their mixing into z> - Since v e = z/ + i> , with (s^) standing for cos( (sin£), 
the survival probability of u e oscillations reads 

P(u e -> u e (t)) = 4 + s 4 ( P(z> -> + 2 c\ s\ Re(u \u (t)), (14) 

where P(vq — > v (t)) is the corresponding "survival probability" for u Q . In terms of the mass 
eigenstates 

N e i(mlt/2E) 

^o(*) = U *, (15) 

n=l Vn 

where we have cut off the sum by N m kvrR, explicitly decoupling the heavy eigenstates. 
This is justified by the arguments given above. For the remaining states, one can show that 
Vn ~ \i^ n ~ I) 2 which follows from the following exact expression for rj 2 : 

r)l = 1 + - csc 2 (27rm„ J R) 2iT 2 m 2 n R 2 + -nsin{2-nm n R) + 2cos(2irm n R) - 2 (16) 

It means that the main contribution to uq comes from the lightest mode, which is expected 
since the bulk zero mode is its main original component. The final survival probability after 
the neutrino traverses a distance L in vacuum can be written down as 



P m {E) = 1-4«E4 sin 2 (^) - 4 4 £ -I, sin= 

„=1 Vn V / k<n VnVk 



(ml - m 2 k )L 



AE 



(17) 



Therefore, the oscillation length is given by L osc = ER\ ~ (E/MeV) x 5 • 10 4 m for Ri = 
0.1 mm. This value for the oscillation length is right in the domain of accessibility of the 
KAMLAND experiment In Fig. 1 we present the survival probability as a function of 
distance for specific values of kvr and hv. Also, we present the probability for the early 
proposal of Ref. || in Fig. 2 for comparision. Notice that the differences in both profiles 
come, basically, from two facts: first, the masses of the light bulk modes have been shifted 
down by |/xo in the present case. Such an effect is absent in the approach followed in the 
former case 0. This is reflected in a (four times) larger oscillation length. Secondly, the 
mixing of v e with bulk states is different in our case than Ref. [0,|J though for large values 
of n they coincide. 
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The averaged probability is now obtained in a stgraightforward manner to be P ee = 
c t l s C' wmc h i s smaller than the two neutrino case with the same mixing angle, although, 



for £ <C 1 it approaches the former result P ee ~ 1 
the modes, but the lowest frequency one, to get 



2£ . Moreover, we may average over all 



p-w = <$ + !«{ +;!<$«? 



1 - 2 sin 2 



1QE 



Hence, the depth of the oscillations is now of the order of ^c 2 c s 2 



(18) 



In conclusion, we have presented a new way to understand small neutrino masses in 
models with large extra dimensions by including the righthanded neutrino among the brane 
particles. In addition to several quantitative differences from earlier works which have only 
then standard model particles in the brane, we find that the string scale in the new class of 
models is necessarily larger. However, we need one extra dimension to be of large size if we 
want to solve the solar neutrino problem as well have other meaningful oscillations between 
the known neutrinos with the bulk neutrinos. 

We have not discussed the matter effect and the implications for solar neutrino puzzle 
nor have we discussed the astrophysical constraints on our scenario. We hope to return to 
these topics subsequently. But we have noted that multiple ripple effect induced by the KK 
modes in the v e survival probability that was noted for the earlier case in Ref. § remains 
in our case too. 

Secondly, we have focussed only on the first generation neutrinos; but there is no obstacle 
in principle to extending the discussion to the other generations. Clearly, as in the case of 
Ref. f|, attempting to solve the atmospheric neutrino puzzle by v^-Vbuik oscillation would 
require that we make the bulk radius R\ smaller. 

Our overall impression is that while it is possible to understand the smallness of the neu- 
trino mass using the bulk neutrinos without invoking the seesaw mechanism and concommi- 
tant high scale physics such as B — L symmetry, a unified picture that clearly incorporates 



all three generations with preferred mixing and mass patterns |14[] is yet to come whereas in 
existing "four-dimensional" models there exist perhaps a surplus of ideas that lead to desir- 
able neutrino mass textures. ^From this point of view, the neutrino physics study in models 
with extra diensions is still in its infancy and whether it grows into adulthood depends much 
on the direction that these extra dimensional models take in the coming years. 
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FIGURES 





/<oL/2E 

FIG. 2. Profile of the survival probability for v e oscillations for the model of first Ref. in [8] 
using £ = 0.1. Notice that the argument is augmented by a factor of four with respect to the case 
in the previous figure. 
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